Impact of an irreversible ß-galactosylceramidase inhibitor on the lipid profile of zebrafish embryos.

Krabbe disease is a sphingolipidosis characterized by the genetic deficiency of the acid hydrolase ß-galactosylceramidase (GALC). Most of the studies concerning the biological role of GALC performed on Krabbe patients and Galc-deficient twitcher mice (an authentic animal model of the disease) indicate that the pathogenesis of this disorder is the consequence of the accumulation of the neurotoxic GALC substrate ß-galactosylsphingosine (psychosine), ignoring the possibility that this enzyme may exert a wider biological impact. Indeed, limited information is available about the effect of GALC downregulation on the cell lipidome in adult and developing organisms. The teleost zebrafish (Danio rerio) has emerged as a useful platform to model human genetic diseases, including sphingolipidoses, and two GALC co-orthologs have been identified in zebrafish (galca and galcb). Here, we investigated the effect of the competitive and irreversible GALC inhibitor ß-galactose-cyclophellitol (GCP) on the lipid profile of zebrafish embryos. Molecular modelling indicates that GCP can be sequestered in the catalytic site of the enzyme and covalently binds human GALC, and the zebrafish Galca and Galcb proteins in a similar manner. Accordingly, GCP inhibits the ß-galactosylceramide hydrolase activity of zebrafish in vitro and in vivo, leading to significant alterations of the lipidome of zebrafish embryos. These results indicate that the lack of GALC activity deeply affects the lipidome during the early stages of embryonic development, and thereby provide insights into the pathogenesis of Krabbe disease.

Myotubularin-related-protein-7 inhibits mutant (G12V) K-RAS by direct interaction.

Inhibition of K-RAS effectors like B-RAF or MEK1/2 is accompanied by treatment resistance in cancer patients via re-activation of PI3K and Wnt signaling. We hypothesized that myotubularin-related-protein-7 (MTMR7), which inhibits PI3K and ERK1/2 signaling downstream of RAS, directly targets RAS and thereby prevents resistance. Using cell and structural biology combined with animal studies, we show that MTMR7 binds and inhibits RAS at cellular membranes. Overexpression of MTMR7 reduced RAS GTPase activities and protein levels, ERK1/2 phosphorylation, c-FOS transcription and cancer cell proliferation in vitro. We located the RAS-inhibitory activity of MTMR7 to its charged coiled coil (CC) region and demonstrate direct interaction with the gastrointestinal cancer-relevant K-RASG12V mutant, favouring its GDP-bound state. In mouse models of gastric and intestinal cancer, a cell-permeable MTMR7-CC mimicry peptide decreased tumour growth, Ki67 proliferation index and ERK1/2 nuclear positivity. Thus, MTMR7 mimicry peptide(s) could provide a novel strategy for targeting mutant K-RAS in cancers.

Comparative analysis of drug-salt-polymer interactions by experiment and molecular simulation improves biopharmaceutical performance.

The propensity of poorly water-soluble drugs to aggregate at supersaturation impedes their bioavailability. Supersaturated amorphous drug-salt-polymer systems provide an emergent approach to this problem. However, the effects of polymers on drug-drug interactions in aqueous phase are largely unexplored and it is unclear how to choose an optimal salt-polymer combination for a particular drug. Here, we describe a comparative experimental and computational characterization of amorphous solid dispersions containing the drug celecoxib, and a polymer, polyvinylpyrrolidone vinyl acetate (PVP-VA) or hydroxypropyl methylcellulose acetate succinate, with or without Na+/K+ salts. Classical models for drug-polymer interactions fail to identify the best drug-salt-polymer combination. In contrast, more stable drug-polymer interaction energies computed from molecular dynamics simulations correlate with prolonged stability of supersaturated amorphous drug-salt-polymer systems, along with better dissolution and pharmacokinetic profiles. The celecoxib-salt-PVP-VA formulations exhibit excellent biopharmaceutical performance, offering the prospect of a low-dosage regimen for this widely used anti-inflammatory, thereby increasing cost-effectiveness, and reducing side-effects.

Molecular Dynamics Simulations of the Human Ecto-5'-Nucleotidase (h-ecto-5'-NT, CD73): Insights into Protein Flexibility and Binding Site Dynamics.

Human ecto-5'-nucleotidase (h-ecto-5'-NT, CD73) is a homodimeric Zn2+-binding metallophosphoesterase that hydrolyzes adenosine 5'-monophosphate (5'-AMP) to adenosine and phosphate. h-Ecto-5'-NT is a key enzyme in purinergic signaling pathways and has been recognized as a promising biological target for several diseases, including cancer and inflammatory, infectious, and autoimmune diseases. Despite its importance as a biological target, little is known about h-ecto-5'-NT dynamics, which poses a considerable challenge to the design of inhibitors of this target enzyme. Here, to explore h-ecto-5'-NT flexibility, all-atom unbiased molecular dynamics (MD) simulations were performed. Remarkable differences in the dynamics of the open (catalytically inactive) and closed (catalytically active) conformations of the apo-h-ecto-5'-NT were observed during the simulations, and the nucleotide analogue inhibitor AMPCP was shown to stabilize the protein structure in the closed conformation. Our results suggest that the large and complex domain motion that enables the h-ecto-5'-NT open/closed conformational switch is slow, and therefore, it could not be completely captured within the time scale of our simulations. Nonetheless, we were able to explore the faster dynamics of the h-ecto-5'-NT substrate binding site, which is mainly located at the C-terminal domain and well conserved among the protein's open and closed conformations. Using the TRAPP ("Transient Pockets in Proteins") approach, we identified transient subpockets close to the substrate binding site. Finally, conformational states of the substrate binding site with higher druggability scores than the crystal structure were identified. In summary, our study provides valuable insights into h-ecto-5'-NT structural flexibility, which can guide the structure-based design of novel h-ecto-5'-NT inhibitors.

PfCRT mutations conferring piperaquine resistance in falciparum malaria shape the kinetics of quinoline drug binding and transport.

The chloroquine resistance transporter (PfCRT) confers resistance to a wide range of quinoline and quinoline-like antimalarial drugs in Plasmodium falciparum, with local drug histories driving its evolution and, hence, the drug transport specificities. For example, the change in prescription practice from chloroquine (CQ) to piperaquine (PPQ) in Southeast Asia has resulted in PfCRT variants that carry an additional mutation, leading to PPQ resistance and, concomitantly, to CQ re-sensitization. How this additional amino acid substitution guides such opposing changes in drug susceptibility is largely unclear. Here, we show by detailed kinetic analyses that both the CQ- and the PPQ-resistance conferring PfCRT variants can bind and transport both drugs. Surprisingly, the kinetic profiles revealed subtle yet significant differences, defining a threshold for in vivo CQ and PPQ resistance. Competition kinetics, together with docking and molecular dynamics simulations, show that the PfCRT variant from the Southeast Asian P. falciparum strain Dd2 can accept simultaneously both CQ and PPQ at distinct but allosterically interacting sites. Furthermore, combining existing mutations associated with PPQ resistance created a PfCRT isoform with unprecedented non-Michaelis-Menten kinetics and superior transport efficiency for both CQ and PPQ. Our study provides additional insights into the organization of the substrate binding cavity of PfCRT and, in addition, reveals perspectives for PfCRT variants with equal transport efficiencies for both PPQ and CQ.

Design, synthesis and biological evaluation of antiparasitic dinitroaniline-ether phospholipid hybrids.

A series of nine novel ether phospholipid-dinitroaniline hybrids were synthesized in an effort to deliver more potent antiparasitic agents with improved safety profile compared to miltefosine. The compounds were evaluated for their in vitro antiparasitic activity against L. infantum, L.donovani, L. amazonensis, L. major and L. tropica promastigotes, L. infantum and L. donovani intracellular amastigotes, Trypanosoma brucei brucei and against different developmental stages of Trypanosoma cruzi. The nature of the oligomethylene spacer between the dinitroaniline moiety and the phosphate group, the length of the side chain substituent on the dinitroaniline and the choline or homocholine head group were found to affect both the activity and toxicity of the hybrids. The early ADMET profile of the derivatives did not reveal major liabilities. Hybrid 3, bearing an 11-carbon oligomethylene spacer, a butyl side chain and a choline head group, was the most potent analogue of the series. It exhibited a broad spectrum antiparasitic profile against the promastigotes of New and Old World Leishmania spp., against intracellular amastigotes of two L. infantum strains and L. donovani, against T. brucei and against T. cruzi Y strain epimastigotes, intracellular amastigotes and trypomastigotes. The early toxicity studies revealed that hybrid 3 showed a safe toxicological profile while its cytotoxicity concentration (CC50) against THP-1 macrophages being >100 µM. Computational analysis of binding sites and docking indicated that the interaction of hybrid 3 with trypanosomatid α-tubulin may contribute to its mechanism of action. Furthermore, compound 3 was found to interfere with the cell cycle in T. cruzi epimastigotes, while ultrastructural studies using SEM and TEM in T. cruzi showed that compound 3 affects cellular processes that result in changes in the Golgi complex, the mitochondria and the parasite's plasma membrane. The snapshot pharmacokinetic studies showed low levels of 3 after 24 h following oral administration of 100 mg/Kg, while, its homocholine congener compound 9 presented a better pharmacokinetic profile.

An engineered variant of MECR reductase reveals indispensability of long-chain acyl-ACPs for mitochondrial respiration.

Mitochondrial fatty acid synthesis (mtFAS) is essential for respiratory function. MtFAS generates the octanoic acid precursor for lipoic acid synthesis, but the role of longer fatty acid products has remained unclear. The structurally well-characterized component of mtFAS, human 2E-enoyl-ACP reductase (MECR) rescues respiratory growth and lipoylation defects of a Saccharomyces cerevisiae Δetr1 strain lacking native mtFAS enoyl reductase. To address the role of longer products of mtFAS, we employed in silico molecular simulations to design a MECR variant with a shortened substrate binding cavity. Our in vitro and in vivo analyses indicate that the MECR G165Q variant allows synthesis of octanoyl groups but not long chain fatty acids, confirming the validity of our computational approach to engineer substrate length specificity. Furthermore, our data imply that restoring lipoylation in mtFAS deficient yeast strains is not sufficient to support respiration and that long chain acyl-ACPs generated by mtFAS are required for mitochondrial function.

Scaling Protein-Water Interactions in the Martini 3 Coarse-Grained Force Field to Simulate Transmembrane Helix Dimers in Different Lipid Environments.

Martini 3, the latest version of the widely used Martini force field for coarse-grained molecular dynamics simulations, is a promising tool to investigate proteins in phospholipid bilayers. However, simulating other lipid environments, such as detergent micelles, presents challenges due to the absence of validated parameters for their constituent molecules. Here, we propose parameters for the micelle-forming surfactant, dodecylphosphocholine (DPC). These result in micelle assembly with aggregation numbers in agreement with the experimental values. However, we identified a lack of hydrophobic interactions between transmembrane helix protein dimers and the tails of DPC molecules, preventing insertion and stabilization of the protein in the micelles. This problem was also observed for protein insertion by self-assembling 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) or dipalmitoylphosphatidylcholine (DPPC) bilayers. We propose the reduction of the nonbonded interactions between protein and water beads by 10% as a simple and effective solution to this problem that enables protein encapsulation in phospholipid micelles and bilayers without altering protein dimerization or the bilayer structure.

Destabilizers of the thymidylate synthase homodimer accelerate its proteasomal degradation and inhibit cancer growth.

Drugs that target human thymidylate synthase (hTS), a dimeric enzyme, are widely used in anticancer therapy. However, treatment with classical substrate-site-directed TS inhibitors induces over-expression of this protein and development of drug resistance. We thus pursued an alternative strategy that led us to the discovery of TS-dimer destabilizers. These compounds bind at the monomer-monomer interface and shift the dimerization equilibrium of both the recombinant and the intracellular protein toward the inactive monomers. A structural, spectroscopic, and kinetic investigation has provided evidence and quantitative information on the effects of the interaction of these small molecules with hTS. Focusing on the best among them, E7, we have shown that it inhibits hTS in cancer cells and accelerates its proteasomal degradation, thus causing a decrease in the enzyme intracellular level. E7 also showed a superior anticancer profile to fluorouracil in a mouse model of human pancreatic and ovarian cancer. Thus, over sixty years after the discovery of the first TS prodrug inhibitor, fluorouracil, E7 breaks the link between TS inhibition and enhanced expression in response, providing a strategy to fight drug-resistant cancers.

MatchTope: A tool to predict the cross reactivity of peptides complexed with Major Histocompatibility Complex I.

The therapeutic targeting of the immune system, for example in vaccinology and cancer treatment, is a challenging task and the subject of active research. Several in silico tools used for predicting immunogenicity are based on the analysis of peptide sequences binding to the Major Histocompatibility Complex (pMHC). However, few of these bioinformatics tools take into account the pMHC three-dimensional structure. Here, we describe a new bioinformatics tool, MatchTope, developed for predicting peptide similarity, which can trigger cross-reactivity events, by computing and analyzing the electrostatic potentials of pMHC complexes. We validated MatchTope by using previously published data from in vitro assays. We thereby demonstrate the strength of MatchTope for similarity prediction between targets derived from several pathogens as well as for indicating possible cross responses between self and tumor peptides. Our results suggest that MatchTope can enhance and speed up future studies in the fields of vaccinology and cancer immunotherapy.

Anti-trypanosomatid structure-based drug design - lessons learned from targeting the folate pathway.

INTRODUCTION: Trypanosomatidic parasitic infections in humans and animals caused by Trypanosoma brucei, Trypanosoma cruzi, and Leishmania species pose a significant health and economic burden in developing countries. There are few effective and accessible treatments for these diseases, and the existing therapies suffer from problems, such as parasite resistance and side effects. Structure-based drug design (SBDD) is one of the strategies that has been applied to discover new compounds targeting trypanosomatid-borne diseases. AREAS COVERED: We review the current literature (mostly over the last 5 years, searched in the PubMed database on 11 November 2021) on the application of structure-based drug design approaches to identify new anti-trypanosomatidic compounds that interfere with a validated target biochemical pathway, the trypanosomatid folate pathway. EXPERT OPINION: The application of structure-based drug design approaches to perturb the trypanosomatid folate pathway has successfully provided many new inhibitors with good selectivity profiles, most of which are natural products or their derivatives or have scaffolds of known drugs. However, the inhibitory effect against the target protein(s) often does not translate to anti-parasitic activity. Further progress is hampered by our incomplete understanding of parasite biology and biochemistry, which is necessary to complement SBDD in a multiparameter optimization approach to discovering selective anti-parasitic drugs.

Identification of antiparasitic drug targets using a multi-omics workflow in the acanthocephalan model.

BACKGROUND: With the expansion of animal production, parasitic helminths are gaining increasing economic importance. However, application of several established deworming agents can harm treated hosts and environment due to their low specificity. Furthermore, the number of parasite strains showing resistance is growing, while hardly any new anthelminthics are being developed. Here, we present a bioinformatics workflow designed to reduce the time and cost in the development of new strategies against parasites. The workflow includes quantitative transcriptomics and proteomics, 3D structure modeling, binding site prediction, and virtual ligand screening. Its use is demonstrated for Acanthocephala (thorny-headed worms) which are an emerging pest in fish aquaculture. We included three acanthocephalans (Pomphorhynchus laevis, Neoechinorhynchus agilis, Neoechinorhynchus buttnerae) from four fish species (common barbel, European eel, thinlip mullet, tambaqui). RESULTS: The workflow led to eleven highly specific candidate targets in acanthocephalans. The candidate targets showed constant and elevated transcript abundances across definitive and accidental hosts, suggestive of constitutive expression and functional importance. Hence, the impairment of the corresponding proteins should enable specific and effective killing of acanthocephalans. Candidate targets were also highly abundant in the acanthocephalan body wall, through which these gutless parasites take up nutrients. Thus, the candidate targets are likely to be accessible to compounds that are orally administered to fish. Virtual ligand screening led to ten compounds, of which five appeared to be especially promising according to ADMET, GHS, and RO5 criteria: tadalafil, pranazepide, piketoprofen, heliomycin, and the nematicide derquantel. CONCLUSIONS: The combination of genomics, transcriptomics, and proteomics led to a broadly applicable procedure for the cost- and time-saving identification of candidate target proteins in parasites. The ligands predicted to bind can now be further evaluated for their suitability in the control of acanthocephalans. The workflow has been deposited at the Galaxy workflow server under the URL tinyurl.com/yx72rda7 .

Diffusion of small molecule drugs is affected by surface interactions and crowder proteins.

Crowded environments are known to affect the diffusion of macromolecules, but their effects on the diffusion of small molecules are largely uncharacterized. We investigate how three protein crowders, bovine serum albumin (BSA), hen egg-white lysozyme, and myoglobin, influence the diffusion rates and interactions of four small molecules: fluorescein, and three drugs, doxorubicin, glycogen synthase kinase-3 inhibitor SB216763, and quinacrine. Using Line-FRAP measurements, Brownian dynamics simulations, and molecular docking, we find that the diffusion rates of the small molecules are highly affected by self-aggregation, interactions with the proteins, and surface adsorption. The diffusion of fluorescein is decreased because of its interactions with the protein crowders and their surface adsorption. Protein crowders increase the diffusion rates of doxorubicin and SB216763 by reducing surface interactions and self-aggregation, respectively. Quinacrine diffusion was not affected by protein crowders. The mechanistic insights gained here may assist in optimization of compounds for higher mobility in complex macromolecular environments.

Combining hypothesis- and data-driven neuroscience modeling in FAIR workflows.

Modeling in neuroscience occurs at the intersection of different points of view and approaches. Typically, hypothesis-driven modeling brings a question into focus so that a model is constructed to investigate a specific hypothesis about how the system works or why certain phenomena are observed. Data-driven modeling, on the other hand, follows a more unbiased approach, with model construction informed by the computationally intensive use of data. At the same time, researchers employ models at different biological scales and at different levels of abstraction. Combining these models while validating them against experimental data increases understanding of the multiscale brain. However, a lack of interoperability, transparency, and reusability of both models and the workflows used to construct them creates barriers for the integration of models representing different biological scales and built using different modeling philosophies. We argue that the same imperatives that drive resources and policy for data - such as the FAIR (Findable, Accessible, Interoperable, Reusable) principles - also support the integration of different modeling approaches. The FAIR principles require that data be shared in formats that are Findable, Accessible, Interoperable, and Reusable. Applying these principles to models and modeling workflows, as well as the data used to constrain and validate them, would allow researchers to find, reuse, question, validate, and extend published models, regardless of whether they are implemented phenomenologically or mechanistically, as a few equations or as a multiscale, hierarchical system. To illustrate these ideas, we use a classical synaptic plasticity model, the Bienenstock-Cooper-Munro rule, as an example due to its long history, different levels of abstraction, and implementation at many scales.

Multitarget, Selective Compound Design Yields Potent Inhibitors of a Kinetoplastid Pteridine Reductase 1.

The optimization of compounds with multiple targets is a difficult multidimensional problem in the drug discovery cycle. Here, we present a systematic, multidisciplinary approach to the development of selective antiparasitic compounds. Computational fragment-based design of novel pteridine derivatives along with iterations of crystallographic structure determination allowed for the derivation of a structure-activity relationship for multitarget inhibition. The approach yielded compounds showing apparent picomolar inhibition of T. brucei pteridine reductase 1 (PTR1), nanomolar inhibition of L. major PTR1, and selective submicromolar inhibition of parasite dihydrofolate reductase (DHFR) versus human DHFR. Moreover, by combining design for polypharmacology with a property-based on-parasite optimization, we found three compounds that exhibited micromolar EC50 values against T. brucei brucei while retaining their target inhibition. Our results provide a basis for the further development of pteridine-based compounds, and we expect our multitarget approach to be generally applicable to the design and optimization of anti-infective agents.

Development and Biological Characterization of a Novel Selective TrkA Agonist with Neuroprotective Properties against Amyloid Toxicity.

Neurotrophins are growth factors that exert important neuroprotective effects by preventing neuronal death and synaptic loss. Nerve Growth Factor (NGF) acts through the activation of its high-affinity, pro-survival TrkA and low-affinity, pro-apoptotic p75NTR receptors. NGF has been shown to slow or prevent neurodegenerative signals in Alzheimer's Disease (AD) progression. However, its low bioavailability and its blood-brain-barrier impermeability limit the use of NGF as a potential therapeutic agent against AD. Based on our previous findings on synthetic dehydroepiandrosterone derivatives, we identified a novel NGF mimetic, named ENT-A013, which selectively activates TrkA and exerts neuroprotective, anti-amyloid-ß actions. We now report the chemical synthesis, in silico modelling, metabolic stability, CYP-mediated reaction phenotyping and biological characterization of ENT-A013 under physiological and neurodegenerative conditions. We show that ENT-A013 selectively activates the TrkA receptor and its downstream kinases Akt and Erk1/2 in PC12 cells, protecting these cells from serum deprivation-induced cell death. Moreover, ENT-A013 promotes survival of primary Dorsal Root Ganglion (DRG) neurons upon NGF withdrawal and protects hippocampal neurons against Amyloid ß-induced apoptosis and synaptic loss. Furthermore, this neurotrophin mimetic partially restores LTP impairment. In conclusion, ENT-A013 represents a promising new lead molecule for developing therapeutics against neurodegenerative disorders, such as Alzheimer's Disease, selectively targeting TrkA-mediated pro-survival signals.

The Non-phosphorylating Glyceraldehyde-3-Phosphate Dehydrogenase GapN Is a Potential New Drug Target in Streptococcus pyogenes.

The strict human pathogen Streptococcus pyogenes causes infections of varying severity, ranging from self-limiting suppurative infections to life-threatening diseases like necrotizing fasciitis or streptococcal toxic shock syndrome. Here, we show that the non-phosphorylating glyceraldehyde-3-phosphate dehydrogenase GapN is an essential enzyme for S. pyogenes. GapN converts glyceraldehyde 3-phosphate into 3-phosphoglycerate coupled to the reduction of NADP to NADPH. The knock-down of gapN by antisense peptide nucleic acids (asPNA) significantly reduces viable bacterial counts of S. pyogenes laboratory and macrolide-resistant clinical strains in vitro. As S. pyogenes lacks the oxidative part of the pentose phosphate pathway, GapN appears to be the major NADPH source for the bacterium. Accordingly, other streptococci that carry a complete pentose phosphate pathway are not prone to asPNA-based gapN knock-down. Determination of the crystal structure of the S. pyogenes GapN apo-enzyme revealed an unusual cis-peptide in proximity to the catalytic binding site. Furthermore, using a structural modeling approach, we correctly predicted competitive inhibition of S. pyogenes GapN by erythrose 4-phosphate, indicating that our structural model can be used for in silico screening of specific GapN inhibitors. In conclusion, the data provided here reveal that GapN is a potential target for antimicrobial substances that selectively kill S. pyogenes and other streptococci that lack the oxidative part of the pentose phosphate pathway.

Cholenic acid derivative UniPR1331 impairs tumor angiogenesis via blockade of VEGF/VEGFR2 in addition to Eph/ephrin.

Angiogenesis, the formation of new blood vessels from preexisting ones, is crucial for tumor growth and metastatization, and is considered a promising therapeutic target. Unfortunately, drugs directed against a specific proangiogenic growth factor or receptor turned out to be of limited benefit for oncology patients, likely due to the high biochemical redundancy of the neovascularization process. In this scenario, multitarget compounds that are able to simultaneously tackle different proangiogenic pathways are eagerly awaited. UniPR1331 is a 3ß-hydroxy-Δ5-cholenic acid derivative, which is already known to inhibit Eph-ephrin interaction. Here, we employed an analysis pipeline consisting of molecular modeling and simulation, surface plasmon resonance spectrometry, biochemical assays, and endothelial cell models to demonstrate that UniPR1331 directly interacts with the vascular endothelial growth factor receptor 2 (VEGFR2) too. The binding of UniPR1331 to VEGFR2 prevents its interaction with the natural ligand vascular endothelial growth factor and subsequent autophosphorylation, signal transduction, and in vitro proangiogenic activation of endothelial cells. In vivo, UniPR1331 inhibits tumor cell-driven angiogenesis in zebrafish. Taken together, these data shed light on the pleiotropic pharmacological effect of UniPR1331, and point to Δ5-cholenic acid as a promising molecular scaffold for the development of multitarget antiangiogenic compounds.